Signal Processing and Control in Nonlinear 
Nanomechanical Systems 



R. L. Badzey 1 , G. Zolfagharkhani 1 , S.-B. Shim 1 , A. Gaidarzhy 1 and 
P. Mohanty 1 

Department of Physics, Boston University, Boston, MA 02215, USA 

1 Introduction 

Bestriding the realms of classical and quantum mechanics, nanomechanical 
structures offer great promise for a huge variety of applications, from com- 
puter memory elements [1] and ultra-fast sensors to quantum computing. In- 
triguing as these possibilities are, there still remain many important hurdles 
to overcome before nanomechanical structures approach anything close to 
their full potential. With their high surface-to-volume ratios and sub-micron 
dimensions, nanomechanical structures are strongly affected by processing ir- 
regularities and susceptible to nonlinear effects. There are several ways of 
dealing with nonlinearity: exceptional fabrication process control in order to 
minimize the onset of nonlinear effects or taking advantage of the interesting 
and oftentimes counterintuitive consequences of nonlinearity. 

A fundamental control mechanism popular for its counter-intuitive ability 
to amplify coherent behavior via the addition of noise, stochastic resonance 
(SR) has waxed and waned in popularity since its inception over 20 years 
ago. Originally, it was postulated [2] as an ad hoc explanation for the pe- 
riodic onset of ice ages over the Earth's climate history. Given a nonlinear 
system with an energy threshold subject to a sub-threshold periodic mod- 
ulation and white noise, there is a certain regime of added noise power in 
which the system oscillates between its states. This oscillation is synchronized 
with the modulation. Although this theory has since fallen out of favor in 
the climate-modeling community, the appeal of the concept encouraged its 
extension to a large variety of systems. These include (but are by no means 
limited to) bistable ring lasers [3] , neurophysiological systems (mechanorecep- 
tors in crayfish [4] and crickets [5]), SQUIDs [6], mechanical systems [7, 8], 
electronic systems [9] such as amplifiers, and geophysical systems [10]. How- 
ever, there have not been any studies demonstrating the effect of stochastic 
resonance in nanomechanical systems. Aside from being simply one more sys- 
tem in which the phenomenon has been demonstrated, nanoscale systems [11] 
are interesting because of their proximity to the realm of quantum mechanics. 
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The combination of stochastic resonance and quantum mechanics has been 
the subject of intense theoretical activities [12, 13, 14, 15] for many years; 
nanomechanical systems present a fertile ground for the study of a broad va- 
riety of novel phenomena in quantum stochastic resonance. Additionally, the 
physical realization of such nonlinear nanomechanical strings offer the possi- 
bility of studying a whole class of phase transition phenomena, particularly 
those modeled by a Landau- Ginzburg quantum string [16, 17] in cosmology. 




2 Nonlinearity in Nanomechanical Structures 

The nanomechanical structure under consideration is a simple suspended 
beam, clamped at both ends and subjected to a transverse driving force. 
Following continuum mechanics, it is well known that such a structure has a 
resonance frequency for the fundamental flexural mode which is given by 



(1) 

where E is the Young's modulus, p is the material density, t is the thickness 
(dimension parallel to the transverse forcing) , and L is the length of the beam. 
Our doubly-clamped beams were fabricated from silicon, with dimensions of 7- 
8 /urn x 300 nm x 200 nm; these yielded resonance frequencies around 25 MHz. 
Clamping losses and electrode mass loading reduced the resonance frequency 
from that which could be assumed from the simple continuum mechanics 
derivation. Fabrication of these structures is a well-established sequence of 
techniques, including PMMA spinning, e-beam exposure, development, met- 
allization, and both dry and wet etching techniques. The details of fabrication 
are given in several locations [1, 18]. The magnetomotive technique has been 
in use for many decades, providing an efficient and low-noise method for ex- 
citing the resonant modes of a suspended structure. A sinusoidal current is 
pushed through a metallic electrode in the presence of a magnetic field. The 
resulting Lorentz force is then: 

'P cos ujt = ~tl COS LOt X ~3 (2) 

It is easy to see that the Lorentz force will be perpendicular to both the field 
and the current; a field perpendicular to the plane of the substrate will produce 
an in-plane transverse vibration. This in turn produces a motional EMF on 
the two clamping electrodes. This voltage is proportional to the velocity of 
the beam: 



V avg (t)=£W^>, (3) 



f dx(t) 
~~dT 

where the factor £ is a proportionality constant given by the mode shape. For 
a fundamental mode, £ = 0.53. In the harmonic approximation, this induced 
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voltage is directly proportional to the displacement of the beam. Under this 
force, the beam will behave as a damped, driven harmonic oscillator: 

mi + "fx + kx = Fcoscdt. (4) 

The amplitude of motion describes a Lorentzian lineshape as a function of 
frequency, centered at the resonance frequency and with a quality factor Q = 
wo/27. In the limit of small dissipation Q = ujo/Aoj. There are two main roads 
into the nonlinear regime: inherent nonlinearities in the suspended bridge, or 
increasing the driving force until the beam response becomes nonlinear. Both 
methods result in the equation of motion described by the Duffing equation, 

m'x + jx + kx ± fc 3 x 3 = F cos tut, (5) 

where the addition of the cubic restoring force term yields a dramatic change 
in the resonance behavior. It is a well-known phenomenon that a Duffing 
oscillator exhibits hysteresis and bistability near the linear resonance fre- 
quency. The evolution of an oscillator from linear response into nonlinear 
behavior is well-studied and covered by a number of excellent texts a sim- 
ple qualitative explanation will suffice. Simply put, as the beam is subjected 
to ever-increasing driving forces, it stretches in response, thereby increasing 
its length. The clamping points being a fixed distance apart, the situation is 
reached whereupon the length of the bridge exceeds that between the clamp- 
ing points. Therefore, when the bridge passes through equilibrium, it feels a 
compressive force from the pads, giving rise to a quartic term in the beam 
potential and therefore a cubic term in the equation of motion. This is the 
well-known Euler instability, which occurs under the influence of any com- 
pressive or tensile force. For a doubly-clamped beam, there is a critical force 
which describes the boundary between the linear and nonlinear regimes [19]: 




Depending on the elasticity of the material, it is possible to switch between 
the linear and nonlinear response regimes at will by simply tuning the driving 
force. Eventually, however, the elastic response gives way to a plastic deforma- 
tion of the bridge and the nonlinear effects become more and more prevalent 
and permanent. When excited nonlincarly, the response of the bridge evolves 
from the simple Lorentzian into curved peak with a dramatic drop as seen in 
the following figure. 

Once the beam is driven into a nonlinear regime it demonstrates hysteresis 
as a function of frequency, as predicted by mechanics. As the frequency is in- 
creased, the amplitude of the response function follows the pseudo-Lorentzian 
lineshape until it reaches the sharp drop as seen in Figure 1 . In reality, the full 
analytic solution of the frequency-response function shows a rather dramatic 
bend to lower frequencies, creating a frequency domain in which the amplitude 
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Fig. 1. Beam response with increasing drive force, showing the transition from the 
linear to nonlinear regime. The labeled curves represent the linear response regime 
of the oscillator in response to different applied forces, from 0.5 pN (curve A) to 3 pN 
(curve D). Beyond this critical force, the beam response is nonlinear and bistable. 




function is multi- valued. Therefore sweeps in one direction in frequency follow 
the upper curve of this region, while opposite sweeps follow the lower curve. 
Eventually, each sweep reaches a frequency in which the system is unstable 
and therefore causes a transition to the other state. This is clearly illustrated 
in Figure 2. 
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Fig. 2. Hysteresis caused by the nonlinear bistable response of a silicon nanome- 
chanical bridge under strong driving. 



As stated previously, the Euler instability occurs under the influence of 
either a compressive or tensile strain on the beam. These two situations lead 
to qualitatively different but dynamically similar situations. In the Duffing 
equation, this is taken into account by the alternative sign on the cubic term, 
with compression yielding a positive sign and tension a negative sign. These 
two situations also present slightly different frequency response spectra as well. 
The compressive case is revealed by a sharp drop on the right-hand side; the 
tensile spring has a drop on the left. For illustrative purposes, Figure 1 shows 
a compressive case, while Figure 2 shows what appears to be a tensile case. 
Regardless of whether there is compression or tension, the oscillator is still 
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Fig. 3. Temperature effects on the dissipation and frequency shift of a nanome- 
chanical silicon bridge, (a). Variations in the resonance peak with temperature. As 
the temperature is increased the peak gets broader and also undergoes shift to lower 
frequencies, (b). Graph of the frequency shift as a function of temperature. The total 
shift is approximately 9.5 kHz, or 4 x 1CT 4 of the initial resonance at 300 mK. (c). 
The dissipation (inverse Q) as a function of temperature. This quantity changes by 
almost a factor of 2 over the temperature sweep. 



bistable and hysteretic. The fact that bistability arises in both nonlinear cases 
is important, because the effect the electrical measurement circuit can have 
on the apparent behavior of the oscillator is worth noting. It is understood 
that a mechanical harmonic oscillator system can be written as an equivalent 
LCR circuit [20] by the same token an electrical system can be constructed 
as a mechanical one. Therefore, when a mechanical system is connected to 
an electrical LCR measurement circuit, the signals measured are due to both 
constituent parts. The electrical resistance, capacitance, and inductance can 
be related to the mechanical spring constant and dissipation 

R~^r, C~k L~\. (7) 

7« k 

These relationships also hold as the oscillator evolves into nonlinear response - 
the electric circuit need not be nonlinear for its signal to combine with that of 
the oscillator. Therefore, the relationship between the effective linear (k) and 
nonlinear (k^) spring constants of the combined electromechanical system can 
be either compressive (same sign) or tensile (opposite sign). Regardless, the 
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physical state of the mechanical structure itself is determined by the physical 
geometry of the system in all of the cases described here, the compressive 
response is the physically realistic one even though the frequency response 
might suggest otherwise. This supposition was borne out by an experimental 
check - changing the cables (to ones with different electrical parameters) re- 
sulted in a flip of the nonlinear response from a left-hand drop to a right-hand 
drop. 
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Fig. 4. Alteration of the hysteresis region by an increase in the temperature from 
300 mK to 4 K. 

In addition, since both tensile and compressive cases result in bistability, 
the distinction from the standpoint of device applications is a largely aca- 
demic one. From the derivation of the frequency response, assuming that the 
solution can be written as a harmonic function, the states are of the vibra- 
tional amplitude. In actuality, the most accurate statement is that they are 
different vibrational velocity states. This is borne out by the response seen by 
the magnetomotive technique, which is at its heart a velocity measurement. 
Only should the harmonic approximation prove valid can it be said that the 
two states are indeed amplitude states. For all of the illustrative abilities of 
the previous analyses, they all neglect the effect of temperature on the oscil- 
latory characteristics of the bridge. At room temperature, slight fluctuations 
will have a negligible effect, and large variations in temperature are neither 
realistic nor easy to implement. Low temperatures, however, are easily at- 
tained through a variety of techniques. We inserted a nanomechanical silicon 
bridge into the cryostat of a 3 He refrigerator and varied the temperature from 
300 mK up to 5 K. 

As is clear from the Figure 3, both the spring constant and the dissipation 
were affected by changes in temperature. The first panel is a series of frequency 
scans taken at a variety of temperatures. Changes in temperature have two 
separate but related effects. First is a shift in the resonance frequency largely 
due to a reduction of the spring constant. The second is a gradual broadening 
of the peak and decrease in the amplitude at resonance, stemming from an in- 
crease in the dissipation factor. However, these two results are not completely 
independent. The frequency shift is also partly responsible for the change in 
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measured Q, making a true separation of the effect of temperature on 7 and k 
difficult. However, to a large degree, measuring Q gives a good approximation 
of 7, and measuring the shift in frequency sheds light on the change in k. 

From 300 mK to 5 K, the frequency shifts by 9.5 kHz; taken as a fraction 
of the initial resonance frequency of approximately 22.0945 MHz, this shift is 
4 x 10~ 4 . This means that the spring constant also changes by approximately 
four parts in 10 4 , as Af ~ Ak. Interestingly, what is also clear from this 
temperature sweep is the effect on the dissipation, measured by tracking the 
changes in the Q of the oscillator. While the frequency shift is negligible, the 
change in dissipation most certainly is not. The dissipation (1/Q) at 300 mK 
has a value of 25 x 10~ 5 ; heating to 5 K increases this value to 45 x 10~ 5 , an 
increase by almost a factor of two. 

It is important to note that the preceding data comes from a nanome- 
chanical oscillator in linear response. In nonlinear response, the notions of 
resonance frequency and Q become very difficult to formulate. However, the 
characteristics of the hysteresis region itself contain information about the 
behavior of 7, k, and k%. These three parameters are even more closely linked 
than in the linear regime. As the concept of resonance frequency is somewhat 
misleading for the nonlinear case, the better measure is the central frequency 
of the nonlinear hysteresis region. By the same token, the concept of Q as the 
ratio between the resonance and the FWHM of the peak loses merit in the 
nonlinear case; a measurement of the change in the width of the hysteresis is 
a better indicator. Again, extracting the exact contribution due to each par- 
ticular parameter is difficult when looking at the temperature alone. When 
the temperature is increased as shown in Figure 4, there are two main results 
a slight shift to lower frequencies and a reduction in the width of the hys- 
teresis. From a low of 300 mK to a high of 4 K, the central frequency of the 
hysteresis region shifts by approximately 1.25 kHz, and the width decreases 
by approximately 500 Hz. It is important to note, however, that transitions 
within the hysteresis region are probabilistic and affected by the ambient tem- 
perature. Many sweeps through the hysteresis region are required in order to 
fully characterize the frequency shift and width change. In general, though, 
this result agrees with earlier experiments [21] that examined the change in 
the nonlinear characteristics with the addition of white noise. 



3 Control by Stochastic Resonance 

One of the more interesting nonlinear effects to come forth in the past few 
decades is the phenomenon of stochastic resonance. Compelling because of 
its ability to draw out coherence from a background of noise, stochastic res- 
onance has suffered a bit in recent years because of a dearth of new and 
interesting experimental observations. Most have been received as simply one 
more example of stochastic resonance in a system, with little impact beyond 
adding to a pile of curiosities. However, the presence of nonlinear behavior in 
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a nanomechanical system [22] opens up an interesting realm of inquiry, a new 
avenue into the exploration of the phenomenon. 

These nanomechanical systems, by virtue of their small sizes, high fre- 
quencies, and low temperatures venture quite close to the regime of quantum 
mechanics. Additionally, these nanomechanical bridges are closely related to a 
system which has become even closer to reality a truly macroscopic, mechan- 
ical quantum harmonic oscillator [23, 24]. One of the very basic and necessary 
conditions for any oscillator to exhibit quantum behavior is that its modal 
energy must be greater than the thermal spreading of the energy states. That 
is, hf > ksT, where ks is the Boltzmann constant and h is Planck's con- 
stant. For an oscillator with a resonance frequency of 1 GHz, this condition 
is achieved when the temperature is below 48 mK. Both of these are experi- 
mentally realizable conditions, and several interesting experiments have been 
put forth which show the deviation of such oscillators from the classical norm. 
Therefore, exhibiting stochastic resonance in these MHz-range nanomechan- 
ical oscillators can lead to exploring the phenomenon in the closely related 
oscillators beyond the quantum limit. At the quantum level, prying coherent 
behavior out of a noisy environment could prove to be a powerful means of 
achieving quantum control. 
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Fig. 5. (a) Schematic of the electrical circuit used to excite a nanomechanical 
silicon oscillator into nonlinear response, modulate switching between the two resul- 
tant bistable states, and introduce electrical broadband noise, (b) Sample switching 
events, showing the effect of reducing modulation. The bottom panel shows the 
synch out from the modulation source. 



To look more universally, these systems closely approximate quantum 
strings described by a Landau-Ginzburg equation. 

d<P(x,t) , , ~. , d 2 <P(x.t) . , 

at at 2 - 

Here, the field variable $ replaces the position variable x, and r\ is the noise 
term. Therefore, these and related nanomechanical structures provide a pow- 
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erful tool for examining fundamental properties of a whole host of systems 
which are governed by similar equations and the resultant phase transition 
effects. The central requirement for the observation of stochastic resonance is 
a nonlinear two-state system subjected to a subthreshold modulation in the 
presence of tunable white noise. With the well-behaved bistable behavior seen 
in these nanomechanical oscillators, they present a very clean system in which 
to study the effect. Also, as the temperature and the driving power affect the 
hysteresis, the beam characteristics can be adjusted so as to study an entire 
parameter space. 

For this experiment, we used two silicon bridges, cooling them in the 3 He 
cryostat and exciting in-plane vibrational modes with the magnetomotive 
technique. At 300 mK bridge 1 had a linear resonance frequency of approxi- 
mately 23.57 MHz (Q = 3700), while bridge 2 had a linear resonance frequency 
of 20.835 MHz (Q = 1000). Both exhibited nonlinear bistable behavior under 
the influence of suitably strong forcing (4.0 dBm for bridge 1, 1.0 dBm for 
bridge 2) . It was found during previous studies of bridge 1 that it was possible 
to excite the oscillator nonlinearly at a single frequency within the bistable 
region. With the addition of a square-wave modulation it was then possible 
to control the state of the oscillator, forcing it back and forth between its two 
states with the modulation frequency fl. 

Under the influence of the driving force and modulation, then, the equation 
of motion for cither of these bridges is a modified Duffing equation of the form: 



This shows the equation with square wave forcing, although in principle any 
modulation term can be used. The response of the bridge was highly dependent 
on the amplitude of the modulation, with low powers leading to a loss of 
switch fidelity but not in the amplitude of the switches. This is consistent 
with expectations, as the voltage difference of the hysteresis is fixed by the 
driving amplitude and temperature the modulation only serves to overcome 
(or fail to overcome) the existing barrier between the two states. Figure 5 
illustrates both the circuit and sample layout and an example of the switching 
behavior as a function of modulation power. Once controllable switching was 
established, the modulation was reduced to sub-threshold, and a broadband 
electrical noise source was introduced into the circuit, and the noise power 
was swept from 71 dBm (~ 79 pW) to 41 dBm (~ 79 nW). At each noise 
power, a time-resolved scan was performed in order to establish whether or 
not the addition of the noise would result in the re-emergence of switching. 
Low noise powers resulted in no response from the bridge it remained in the 
initially prepared state. As the noise was increased, however, first sporadic, 
and then more synchronized switch events were seen. Each switch event was 
synchronized with the modulation, and occupied exactly one-half period of 



mx + -fx + kx ± fc 3 x 3 = Fd cos ui^t + / m (4?) 
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the modulation frequency of 0.05 Hz. If a switch was skipped, no new switch 
would appear until the modulation passed through at least one-half period 
and again entered the positive-going part of the signal. The initial increase in 
switching was dramatic and quite rapid, establishing full switching within only 
a few fractions of a dBm increase in noise power. The signature of stochastic 
resonance is an increase in the signal-to-noise ratio (SNR) as a function of 
applied noise power. SNR is defined as: 

here, S(f2) is the height of the peak in the power spectral density at the 
modulation frequency, and N(f2) is the background noise. With the SNR 
data accumulated, the following graph was compiled, showing a quite dramatic 
increase in the SNR with increasing noise power. An interesting feature of this 
data is the fact that the SNR does not exhibit the typical sharp peak as a 
function of the noise, but is instead rather broad, extending for several dBm 
before declining sharply down to near zero again. It should be noted that 
our system was distinct for several reasons. One, it was subjected to a square 
wave modulation as opposed to the canonical sine wave. The effect of a square 
wave on stochastic resonance has been the subject of considerable theoretical 
interest [26, 25]. Secondly, the electronic noise source was of broad but still 
finite bandwidth - 15 MHz, to be exact. Canonical theory quite explicitly 
states that the noise introduced is Gaussian and broadband. 
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Fig. 6. SNR as a function of applied noise power, (a) Three representative time 
series, showing the reemergence of switching as the noise power is increased. The 
applied noise increases from the uppermost panel downward, (b) Taking a power 
spectrum of each time series yields the SNR for each noise power these were com- 
bined to create a graph of the SNR versus the applied noise power. The signature 
increase of stochastic resonance is clear, even though the overall shape does not 
follow the expected canonical model. 
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Additionally, because of the very clean nature of the system switching, 
the barrier between the two states was quite strong. There is a considerable 
amount of rigidity to the system - it is difficult to induce switching between 
the two and the switching, once begun, is robust. By the same token, dynamic 
changes to the conditions of the oscillator were not felt instantaneously. It is 
possible that a longer time series at each noise power would reveal subtle 
differences in the switching behavior that would lead to changes in the SNR. 
Finally, the nature of the oscillator also lends itself unusual switching charac- 
teristics once the noise has increased above the effective power for stochastic 
resonance. In the canonical formulation and most experiments in stochastic 
resonance, the decrease in the SNR at higher noise powers is due to a swamp- 
ing of the switching signal by a large number of incoherent switches. However, 
this is not the case here the oscillator simply stops switching and instead re- 
mains in the upper state. As this is a symmetric system, there is no preference 
for the eventual settling in the upper state - subsequent runs would leave the 
system in the lower state with just as much probability. And lastly, the effect 
of the fact that the noise source bandwidth did not intersect or overlap with 
the system resonance frequency should be explored in more detail, as it could 
contribute to the rigidity. This matter of system rigidity is certainly an open 
question, and would benefit from additional experimentation and theoretical 
consideration. 

As temperature has already proven itself to be a powerful noise source [18] 
in the behavior of nanomechanical oscillators, a natural check was to see if 
stochastic resonance occurred with the source of noise being temperature, not 
external and electrical. 




Fig. 7. SNR as a function of ambient temperature, (a) A sequence of time series 
for this asymmetric bridge shows that even when the coherence is regained, it is not 
as clean as for the first bridge. Also, the first panel, which shows the sharp spikes 
and immediate decays, demonstrates the asymmetric character of this bridge, (b) 
The SNR graph demonstrates a strong resonance. The data is simply line-connected 
there is no fit implied, and is affected by the coarseness of the sweep. 
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The second oscillator was disconnected from the noise source and brought 
into nonlinear response. The modulation required for switching was excep- 
tionally strong (19 dBm), to the extent that the electrical signal appeared in 
the time series on top of the normally-static beam response voltage. Again 
reducing the modulation to below threshold (16.5 dBm), the temperature was 
increased from 300 mK to 4 K and the SNR was extracted from the time 
series of the beam response. It is clear from Figure 7, there is another peak in 
the SNR as the temperature increases beyond 2 K. Although it is tempting 
to apply the same analysis to the SNR response of the second bridge as the 
first, the results should be scrutinized carefully. It is obvious from the time 
series and the power spectra that this bridge had a non-negligible asymme- 
try. The presence of asymmetry has been broached in a number of theoretical 
investigations [27], and it can be used to understand this system, as well. 

Because the temperature is changing, its effect on the dissipation and 
spring constant of the oscillator should not be overlooked. As shown in the 
first section, the dissipation is quite strongly changed as a function of the 
temperature. Therefore it is necessary to look at the equation of motion with 
new eyes. Typically, the noise considered in the formulation of the stochas- 
tic resonance problem is purely additive. However, when the dissipation is a 
function of temperature, the damping coefficient 7 can change the dynamics 
of the situation. In this case, it would be best to consider the noise as being 
partially additive (from the direct effect of heating brought on by the increase 
in temperature) and partially multiplicative (from the changes to dissipation) . 
Over the temperature range where the stochastic resonance is seen (2 K to 

3 K), the dissipation (based on the linear data) changes by approximately 10 
per cent - certainly it is an effect worthy of further investigation. The effect 
on the spring constant, however, is not large. 

In conclusion, we have introduced stochastic resonance into an entirely new 
class of systems, ones which by virtue of their small size and high frequencies 
touch upon an area where this and other nonlinear phenomena may have 
a deep and lasting impact. Nanoscale systems, and nanomechanical systems 
in particular, are much more than another example of stochastic resonance. 
Their highly controllable nature makes them perfect candidates for the study 
of a host of exciting phenomena, to an extent rarely seen before. 
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